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Abstract 

Equations of geodesic deviation for the 3-dimcnsional and 4-dimensional 
Riemann spaces are discussed. Availability of wide classes of exact so- 
lutions of such equations, due to recent results for the matrix Schrodinger 
equation, is demonstrated. Particular classes of exact solutions for the 
geodesic deviation equation as well as for the Raychaudhuri and gen- 
eralized Raychaudhuri equation are presented. Solutions of geodesic 
deviation equation for the Schwarzshild and Kasner metrics are found. 



1 Introduction 

It is well-known that the equation of geodesic deviation is an important 
equation in general relativity. It relates the relative acceleration between 
two test particles to certain components of the Riemann-curvature tensor. 
A knowledge of geodesic deviations is also needed to evaluate or simplify the 
first and higher derivatives of two-point geometrical quantities( the world 
function or geodetic interval), the parallel propagator and others geometric 
quantities related to these [1-3]. 

The Raychaudhuri equation, which is closely connected with the equa- 
tion of geodesic deviation plays a key role in an analysis of the focusing effects 
of gravity [3,4]. Evolution of deformations of a relativistic membranes and 
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the focusing of families of surfaces are described by the generalisations of 
the Raychaudhuri equation [5,6]. 

Equation of geodesic deviation is also an important in mathematics, in 
the theory of Riemann spaces. For surfaces it coincides, in essence, with 
the Gauss equation in geodesic coordinates while in the three and higher- 
dimensional Riemann spaces it carries an essential information about these 
spaces. Solutions of the geodesic deviations equation (Jacoby fields) and 
their properties(e.g., an existence of conjugate points) are related to various 
important characteristics of Riemann spaces [7-8]. 

The equations of geodesic deviation are the system of second order linear 
equations. Their relevance to various problems in physics and mathematics 
has became clear long time ago. But, suprisingly, still not many exact 
solutions of these equations have been found. They basically correspond to 
the particular case of reduction to a single second order equation [4,6,9]. 

In this paper we will present a wide variety of exact soliutions of the 
equation of geodesic deviations. These solutions became available due to 
the recently discovered method of the inverse spectral transform [10,11]. 
We will consider equations of geodesic deviation in the two, three and four 
dimensional Riemann spaces as well as the Raychaudhuri equation and its 
generalisations. 

The case of the diagonalized separated system of the geodesic deviation 
equations which correspond to the diagonal metric is considered in details. 
The nondiagonalizable case is studied with the use of the theory of the 
matrix Schrodinger equation. General formulae and particular solutions are 
presented. 

Specification of the geodesic deviations equation via nonlinear integrable 
(soliton) equations is discussed. Solutions of the geodesic deviations equa- 
tion in general relativity for the Schwarzschild and Kasner metrics are pre- 
sented. 

2 Equation of geodesic deviations 

Here we will present some known formulas, which will be used in what 
follows. In the Riemann space with the metric 

ds 2 = gijdx i dx^ (1) 
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the geodesies are defined via the equation 

d?x l i dx l c 
~d^ +Tkj lh~ds 



d 2 x l z dx % dxi , . 

n,-r--r- = o (2) 



where Y l k - are the Christoffel symbols of metric. 

Let rf be the components of deviation vector between two infitesimally 
nearby geodesic lines. Then these components satisfy to the Jacoby equation 
[1] 

jViVVrf) = -w'i&^V (3) 

where v l are the components of tangent vector along geodesic line 7, symbol 
V, means the covariant derivation along the vector field v l with respect to 
the metric (1) and Rj kl are the curvature tensor of the metric (1). 

In a special system of coordinates were axis x % is a geodesic line equations 
(3) have the following form [1-3] 



d 2 rf , , 



+ RLri l = 0. (4) 



This system of equations is a main object of our consideration. 

The system of equations (4) has a form of the matrix Schrodinger oper- 
ator: 

[-1^ + U(x)]9 = \9 (5) 

where Uij(x) = (—R J iH + XSji), x l = x and * is the fundamental matrix- 
solution. 

In the one-dimensional case the equation of geodesic deviation is equiv- 
alent to the Gauss equation 

fl+K(x,y)* = (6) 

which connects the curvature K(x, y) of the surface with the metric 

ds 2 = dx 2 + ^ 2 {x,y)dy 2 . (7) 

The Raychaudhuri equation also is, in fact, the one-dimensional Schrodinger 
equation [4,3]. Indeed, introducing the variable = JMogdetC/, one gets 
from equation (4) the following equation 

d@ 1 ^ 9 9 „ 
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e = 3-~ x 



where a 2 = a ab a a i), u) 2 = u ab oj a b and R are certain quantities depending on 
the metric. The change of variable 

leads to the equation 

V xx = ^LU 2 -a 2 - R)V. 
In particular, for the isotropic cosmological model with the Einstein equation 

-^T» = R»- l -R5t + k5t 

one has 

|^ = i[A-47rp-^]* (8) 

where p is a density of matter and is a certain geometric characteristic. 

The Raychaudhuri equation is the principal equation for an analysis of 
the singularity behaviour in general relativity. Few particular solutions of 
this equation have been found untill now. 

Wide classes of solutions of equation(6) have been presented in the paper 
[12] . Since for surfaces r\ = ^(x,y)dy then the results of [12] provide wide 
classes of the Jacoby field for the surfaces. 

3 Geodesic deviation equation for three-dimensional 
space 

In the case of three-dimensional space the use of the coordinate transforma- 
tions allows us to reduce a general metrics to the form 

ds 2 = dx 2 + A(x, y, z)dy 2 + 2B(x, y, z)dydz + C(x, y, z)dz 2 . (9) 

For such a metric the matrices of the ChristoffeFs symbols are 



n CA X —BB X CB X —BC X 
u 2A 2A 



V o 



AB X —BA X AC X ~BB X 



2A 



2A 



4 



r 2 = 



Ax 
2 



2 



CA X ~BB X C Ay+BAz~2BB y CA z ~BC y 
2A 2A 2A 

AB X -BA X 2ABy-BA y -AA z AC y —BA z 



2A 



2A 



.2 

2A / 



B x 
2 



2 



CB X ~BC X CA z -BC y 2CB z -CC y -BC z 
2A 2A 2A 



AC X -BB X AC y —BA z AC z +BC y -2BB z 
\ 2A 2A 2A 



where A = AC - B 2 . 

The equations of geodesic deviations are of the form 



dx 2 
dx 2 



+ R 2 2lV 2 + R 2 zl rf = 0, 
+ Rl 2lV 2 + R 3 13lV 3 = 



(10) 



where the corresponding components of the curvature tensor are 

1 d . CA X - BB X . (CA X - BB X ) 2 (CB X - BC X ) {AB X - BA X ) 



R2l21 ~ 



^131 - - 



A ' 4A 2 4A 2 

1 d CB X — BC X (CB X — BC X )A X 



2dx 



A 



-)- 



4A 2 



3 1 8 ( AB X -BA X ^ (AB X - BA X )A X 

Kl21 ~ 2d^ { A ] IK 2 ' 

3 Id AC X - BB X {AC X - BB X ) 2 (AB X - BA X )(CB X - BC X ) 

131 2dx { A ' 4A 2 4A 2 

The system of equations (10) can be rewriten in form of the matrix 
Schrodinger operator 



d 2 V 2 
dx 2 
d 2 r\ z 
dx 2 



+ (-R 2 21 + \ 2 )r ] 2 + (-R 2 31 W = \ 2 r ] 2 , 
+ {-R\2i)ri 2 + (-i2? 3 i + A V = \W- 



(11) 
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4 Some examples of solutions of geodesic devia- 
tion equation 

The simplest case of the deviation equations (7) corresponds to the diagonal 
metrics, i.e. to the case B = 0. In this case we have the system 

dx 2 ^ dx 2 ^ ^ ^ 

dW d 2 B 3 n 
— H rj = 

dx 2 dx 2 

and its solutions correspond to the geodesic deviations for the 3-dimensional 
orthogonal metrics (9). We can construct wide classes of solutions of each of 
equations (12) using the results for the one-dimensional Schrddinger equa- 
tion. In this case the components of deviations vector rj 2 and if are com- 
pletely independent. 

One can relates them by some special conditions. One possibility is to 
consider the system of equations (12) as related by the Darboux transfor- 
mations, i.e. 

d 2 B d 2 A „ d 2 



dx 2 dx 2 dx 2 

where rj 2 is some solution of the first equation (12). The set of metrics 
related by the conditions of such type remains always undetermined and for 
its further specification it is necessary to fix the dependence on the variables 
y and z. We can fix this dependence with the help of the some equations or 
to consider the metrics as induced from the immersion in some Euclidean 
space or another type of spaces. 

As example we can consider the class of 3-dimensional orthogonal metrics 

ds 2 = A 2 (x, y, z)dx 2 + B 2 (x, y, z)dy 2 + C 2 (x, y, z)dz 2 

and will require that the functions A, B, C obey the Darboux equations 
[13,14] 

. Cy . B z 

A zy = "^tAs + ~q y> 



B Z x — ~prB z -\ -r-B x , 
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r - ^Lr a- —r 
y ~ A B 

In the case A = 1 this system is reduced to the relations 

B z = U(y,z)C, C y = V(y,z)B, 

which are equivalent to the Laplace-type equation for the functions A and 
B 

B zy - ^-B z - UVB = 0, C zy - ^rC y - VUC = 0. 

These systems of equations are well known in theory of integrable equations 
and one can use different methods to construct their exact solutions. 

In any case using the N-times iterated Darboux transformation and cer- 
tain addition conditions one can get the discrete set of metrics 

ds 2 n = dx 2 + B n (x,y,z)dy 2 + C n (x,y, z)dz 2 , 

which may be of interest from the various points of view. 

For analysis of nondiagonal metrics (9) we can use the results on the 
matrix Schrodinger equation and corresponding nonlinear integrable differ- 
ential equations. 

The simplest way consists in the use of the relations between the AKNS- 
system [10,11] 

dip! 



dx 
dip 2 



- + iAV>i = q(x,y,z)i/j 2 , 
i\ip 2 = r(x,y,z)i/)i 



dx 

and the corresponding Schrodinger-like system of equations [15]: 



1 \ cP_ + / rq q x 
1/ dx 2 I r x rq 



t 1 = A M t 1 

V2 \ V>2 



Comparing these equations with equations (8), one obtains the following 
conditions on curvature tensor 

A 2 - R\ 2l = rq, A 2 - R 3 131 = rq, 

R-131 = ~Qx, Rl21 = —fx- 

These conditions imply certain relations between components of metric. 
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In fact, from the condition 



-^121 — ^131 



we have 

CA X - AC X = S(y, z)\/AC-B 2 . 
Other two conditions lead to the relations 

AB X - BA X = 2\/ AC - B 2 {^J r x \J AC - B 2 dx + a(y, z)j , 
CB X - BC X = 2V AC - B 2 (^j q x VAC-B 2 dx + 0(y, z 



where a(y,z), (3(y,z) and 5(y,z) are arbitrary functions. 

So, we have some relations between the functions and we can hope that 
the solutions of this system exist. As result, we can obtain the two para- 
metrical family of solutions of the matrix Schrodinger equation and then we 
can fix dependence on variables y and z by nonlinear integrable equation 
which have the form of the Lax-equations 

L y = [L, A] or L z = [L, A] 

where 

d 2 

L = -^2 +U(x,y,z) 
and A is some matrix differential operator. 

So this approach to the solutions of geodesic deviation equations allows 
us to calculate the examples of 3-dimensional metrics which could be useful 
in the various problems (for example, in general relativity). 



5 Solutions of the matrix Schrodinger equation 

and geodesic deviation equation in the 4-dimensional 
space 

We see that the construction of solutions of the matrix Schrodinger equation 
is a key step in the study of the geodesic deviation equation. 

The matrix Schrodinger equation has been studied recently within the 
framework of the inverse spectral transform method [15,16]. We will present 
here some results from [16]. 
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Solution of equation (5) is given by the formulas 

roo 

* = exp(ifcx) + / K(x,x') ex.p(ikx')dx' 

Jx 

where k 2 = X and 



U(x) = -2-^-K(x,x). 

ax 



Here the matrix function K(x, x') is the solution of linear integral equation 

roo 

K(x,x') + M(x + x') + dx"K(x,x")M(x' + x") = 0, x < x' 

Jx 



with the kernel 



N l r + oo 

Mix) = C n exp(—k n x) H / dXR(X) exp(iAx) 

where C n and P(A) are arbitrary matrices. The simplest solution of this inte- 
gral equation corresponds to the case R = and N = 1. The corresponding 
potential £7 and wave function * are 

Q(x) = -2p 2 cosh 2 [p(x - £)]P, (13) 

* = (ifc — ptanh[p(x — £)]) exp(ifc(x — £)).P 
where P is projector, i.e. the matrix which satisfies the condition 

P 2 = P. 

More complicated solution corresponds to N = 2 and R = 0. It is of the 
form 

U(x,y,z) = -W x , 
W{x,y,z) = -2{ Pl + P2 )[l - pt 1 t 2 }- 1 [t 1 P 1 +t 2 P2-t 1 T2{Pi,P2}}, 

Tk = 7 ^ M ^ ~ tanh \Pk(x ~ &)]) 

(Pl+P2jJ 

and the function p is defined from the relations 

pP 1 = P1P2P1, 
pP 2 = P 2 PiP 2 - 
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All solutions which correspond to the case R = can be written in 
a closed form [15,16]. These exact solutions could be used to specify the 
components of the curvature tensor. For instance, for the solution (13) one 
has the following relations between components of curvature tensor 

n i3i n m — - ft i2i- ft i2i- 
In the case of 4-dimensional space with the geodesic coordinate system 

ds 2 = dt 2 + g ab dx a dx b 
the geodesic deviations equation has the form 
d 2 

-^r] 1 + R l mo r, 1 + i4 r/ 2 + R\ m if = 0, 
d 2 

^v 2 + RlioV 1 + R 2 2oV 2 + R 2 3oV 3 = o, 
d 2 

-^if + ^oio^ 1 + RmoV 2 + Rm^ = 

where Rj kl are the components of curvature tensor. In this case the condi- 
tions on Rj kl to be projector have more complicated form. 

Using the generalizations of Wronskian relations between the various 
solutions of the matrix Schrodinger equation one can get the nonlinear 
relations for the potentials: Q(x, y, z) and Q(x, y + Ay, z + Az) at the fixed 
value of variable x such that the corresponding matrix coefficients R(X, y, z, ) 
and C(y,z,x) = SS T will satisfy linear equations. 

In a such a way the nonlinear "boomeron" integrable equations of the 
form 

Q y (x, y, z) = y, z),Q z (x, y, z)] 

have been constructed [16]. Here $ is a certain function. Their solutions can 
be used for the construction of the solutions of geodesic deviation equation. 

So we see that the soliton structures and equations may be useful in the 
study of the geodesic deviations equations and in the general relativity. 

6 Some solutions of geodesic deviation equation in 
general relativity 

Let us consider the case of 4-dimensional Riemann space with the metric 
[17] 

ds 2 = dt 2 - A 2 (t)dx 2 - B 2 (t)dy 2 - C 2 (t)dz 2 . 
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The corresponding geodesic deviation equations are of the form 

d 2 i A tt i n 

dfi 71 -IT" =0 ' 

2 Btt 2 n 
^ =0 ' 

d 2 3 C« 3 n 

If the metric (25) satisfies the Einstein equations, then: 
A = t Pl , B = t P2 , C = t p:i 

where 

Pi+P2+P3 = l, P 2 +P2+P 2 = l - 
It is the Kasner metric. The variables pi admit the parametrisation: 

-s s(s + l) 1 + s 



Pi = ~, — ; ; — o P2 = ^— : — n P3 



1 + s + s 2 ^ 1 + s + s 2 ^ J 1 + s + s 2 ' 
Hence, the corresponding equations are not independent, namely: 

dt 2 " 1 {l + s + s 2 )t 2V ' 

d 2 2 S ( S + 1) 2 

^ + (l + s + s 2) t 2 T l =0 ' 

(i 2 3 (1 + S)s 2 3 _ 

d^ 77 + (l + s + s 2 )^ 
Solutions of the equation of the type 

y" = caT 2 j/ 

have the form: 

i I „ 

C\X 2 + C2X2 , 



U— Cx^/x + C 2 \/xlogX 

ciy / xcos(s log x) + C2y / xsin(slogx), 

where 2s = + 1| and for the first solution 4c + 1 > 0, for the second 
solution 4c + 1 = and for the third 4c + 1 < 0. These solutions provide 
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us the solutions and properties of the geodesic deviations for the Kasner 
metric. 

To study the properies of geodesic deviation for the Schwarzschild solu- 
tion of Einstein equation it is convenient to use the metric in the Lemeitre- 
form [17] 



ds 



2 _ 



dt 2 - 



dR 2 




2/3 


(«"*)_ 





§(12 - t)} 4/3 r 2 g /3 {d9 2 + sin 2 9d(j) 2 ). (14) 



The corresponding equations for geodesic deviations are 

d 2 1 4 
W ~ — 



dt 2 
dt 2 



v 2 + 



v 3 + 



9(12- 


t) 2 


2 




9(12- 


t) 2 


2 




9(12- 


t) 2 



ri l =0, 
V 2 =0, 
r? 3 = 0. 



The solutions of these equations have the following simple form 

rj 1 = Cl (R-t)^ 3 + c 2 (R-t)- 1 / 3 , 

~ 2 -cz(R-t) l l 2 + c A (R-t) l l\ 
c^R-t) l l 2 + c Q (R-t) l l\ 



V 



These solutions can be useful for the study the properties of Schwarzshild 
metric. 



7 Solution of the Raychaudhuri equation 

The Raychaudhuri equation (8) is the one-dimensional scalar Schrodinger 
equation with the potential 

u = Airp + (f> 2 

and the cosmological constant A plays a role of energy (A = ^) So, any ex- 
act solution of the Schrodinger equation provides us the solution of the Ray- 
chaudhuri equation. Wide class of them is given by the , so-called, Bargmann 
potentials and their eigenfunctions 

U(x,y) = 2(logdetA) ;E;E 
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and 

, T , of r -\ i , ^ detA n exp[-(a n (y) +i\ )x] 
w = Re(exp\—i\ x\ + > — ) 

where y denotes all independent variables except x, A is the N x N matrix 
with elements 

A nk = 5 nk H exp[-(a„(y) + a k (y))x]. (15) 

The matrix elements of the matrix A n are given by (15) with the sub- 
stitution of the last column by the column 

-Pn(y) exp(-a n (y)x) 

(n=l,2,3,...N). 

In the simplest case N = 1 one has: 

v = , 2 ° 2to) ± 

cosh [a{y)x — j(y)] 3 

* = ite( iAo + " iM expHAox]) 

zA + a(y) 

where 

, x i , /%) 



2 °2a(y) 

and a(y), /3(y) are arbitrary functions. These solutions have very special 
properties: they are transparent for all energies, the corresponding function 
^ have in general N — 1 zeros. 

The dependence of the functions a, (3 can be fixed by the requirement 
that U and obey additional equations. Theory of soliton equations pro- 
vides us such equations. They are the famous, Korteweg-de Vries (KdV) 
equation and its higher partners. The KdV equation has the form 

U y + 6UU X + U xxx = 

while ^ obeys the equation 

*» + 4® xxx - QU^ X -W x ^ = 0. 

The formulae (20) with a n = const , (5 n = /3„o exp(8a^y) give us the multi- 
soliton solutions of the KdV equation. These solutions have a number of 
interesting properties. 

The Raychaudhuri equation inherits all properties of these exact solu- 
tions. We will discuss the properties of such soliton like solutions within the 
general relativity elsewhere. 
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8 Solutions of the generalized Raychaudhuri equa- 
tion 



A generalisation of Raychaudhuri equation in the case of the two-dimensional 
time-like surfaces embedded in a four-dimensional background has been de- 
rived in [5]. It has the form 

A7+^ fl 7 +(M 2 )i = 0, 

where A = V a V a , V a is the world-sheet covariant derivative and d a ~f = 
a . The quantity {M 2 )\ is connected with the geometric characteristic of 
background and geometry of membrane embedded in this background. 
In a simplest case this equation looks like [6] 

-^ + ^ + ^ (a ' r)(M2)Ka ' r)Xl/ = 

where Q 2 is the conformal factor of the induced metric written in isothermal 
coordinates. 

It is is the second-order linear, hyperbolic, partial differential equation 
for \I> describing the deformation of the surfaces with respect to parameter r. 
It coincides with deviations equation in the case when ^ does not dependent 
on r. 

We present this equation in form 

- (d 2 -d^ + U(T,a)^ = (16) 

where 

U = n 2 (M 2 )l + E. 

Particular solution of this equation has been obtained in [7] by separation 
of variables. A wide class of solutions of equation (16) can be constructed 
by the inverse spectral transform method. The simplest from them are of 
the form [18,19]. For E / 



*(t, a) 



a(3 cosh 2 [(a - (3)(r + a) + + 1 \/2 

cosh[(a + (3){a + r + + S)}/2 

cosh[(a-/?)(a + r + ^^ + 7 )]72 
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where a, j3, 7, 5 are arbitrary parameters. 
For the case i? = 0we have the solution 

{T,(J) cosh 2 [a(r + a) + 0(r - a)] 

and 

q,( T a \ = ^ cosh[g(r + a) - (3{t - a) + 5] 
cosh[a(r + a) + (3{t - a) + 7] ' 
A very simple solution with E = is of the form 

6(T + CT -a(r-a)) 2 -2/3 2 
l ' J (r + a-a(r- CT ) 2 + /3 2 ) 2 

and 

*(t,(t) 



(r + a - a(r - a)) 2 + (3 2 

with arbitrary parameters a, (3. Such exact solutions of could be important 
for an analysis of dynamics of surfaces. 
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